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1 Introduction 

When a fluid is in contact with another fluid, or with a gas, a portion of the 
total free energy of the system is proportional to the area of the surface of 
contact, and to a coefficient, the surface tension, which is specific for each 
pair of substances. Equilibrium will accordingly be obtained when the free 
energy of the surfaces in contact is a minimum. 

Suppose that we have a drop of some fluid, b, over a flat substrate, w, 
while both are exposed to air, a. We have then three different surfaces 
of contact, and the total free energy of the system consists of three parts, 
associated to these three surfaces. A drop of fluid b, will exist provided its 
own two surface tensions exceed the surface tension between the substrate w 
and the air, i.e., provided that 

T wb _|_ T ba > T wa 



If equality is attained, then a film of fluid b is formed, a situation which is 
known as perfect, or complete wetting. 

When one of the substances involved is anisotropic, such as a crystal, the 
contribution to the total free energy of each element of area depends on its 
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orientation. The minimum surface free energy for a given volume determines 
then the ideal form of the crystal in equilibrium. 

It is only in recent times that equilibrium crystals have been produced 
in the laboratory, first, in negative crystals (vapor bubbles) of organic sub- 
stances. Most crystals grow under non-equilibrium conditions and is a sub- 
sequent relaxation of the macroscopic crystal that restores the equilibrium. 
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Figure 1. Partial and complete wetting 

An interesting phenomenon that can be observed on these crystals is the 
roughening transition, characterized by the disappearance of the facets of a 
given orientation, when the temperature attains a certain particular value. 
The best observations have been made on helium crystals, in equilibrium 
with superfluid helium, since then the transport of matter and heat is ex- 
tremely fast. Crystals grow to size of 1-5 mm and relaxation times vary from 
milliseconds to minutes. Roughening transitions for three different types of 
facets have been observed (see, for instance, Wolf et al, 1983). 

These are some classical examples among a variety of interesting phe- 
nomena connected with the behavior of the interface between two phases in 
a physical system. The study of the nature and properties of the interfaces, 
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at least for some simple systems in statistical mechanics, is also an inter- 
esting subject of mathematical physics. Some aspects of this study will be 
discussed in the present article. 

We assume that the interatomic forces can be modelled by a lattice gas, 
and consider, as a simple example, the ferromagnetic Ising model. In a 
typical two-phase equilibrium state there is a dense component, which can 
be interpreted as a solid or liquid phase, and a dilute phase, which can be 
interpreted as the vapor phase. Considering certain particular cases of such 
situations, we first introduce a precise definition of the surface tension and 
then proceed on the mathematical analysis of some preliminary properties of 
the corresponding interfaces. The next topic concerns the wetting properties 
of the system, and the final section is devoted to the associated equilibrium 
crystal. 

2 Pure phases and surface tension 

The Ising model is defined on the cubic lattice C = Z 3 , with configuration 
space f2 = { — 1, 1} £ . If a G fi, the value a(i) = — 1 or 1 is the spin at the 
site % = (ii, «2, ^3) G £, and corresponds to an empty or an occupied site in 
the lattice gas version of the model. The system is first considered in a finite 
box Ac£, with fixed values of the spins outside. 

In order to simplify the exposition we shall mainly consider the three- 
dimensional Ising model, though some of the results to be discussed hold in 
any dimension d > 2. We shall also, sometimes, refer to the two-dimensional 
model, being then understood that the definitions have been adapted in the 
obvious way. We assume that the box A is a parallelepiped, centered at the 
origin of £, of sides L 1; L 2 , L 3 , parallel to the axes. 

A configuration of spins on A (a(i),i G A), denoted <ta, has an energy 
defined by the hamiltonian 

H A (a A \a) = -J ]T a(i)a(j), (1) 

<i,j)nA^0 

where J is a positive constant (ferromagnetic or attractive interaction). The 
sum runs over all nearest neighbor pairs C £, such that at least one 

of the sites belongs to A, and one takes a(i) = a(i) when % £ A, the config- 
uration o G f2 being the given boundary condition. The probability of the 
configuration <ta, at the inverse temperature (3 = 1/kT, is given by the Gibbs 
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measure 



/i A (a A | a) = Z 5 (A)- 1 exp ( - 0H A (a A \ a)) 
where Z a (A) is the partition function, 

^(A) = 5>xp(-/3# A (a A |a)). 



(2) 



(3) 



Local properties at equilibrium can be described by the correlation functions 
between the spins on finite sets of sites, 



The measures (j2J) determine (by the DLR equations) the set of Gibbs 
states of the infinite system, as measures on the set Q of all configurations. 
If a Gibbs state happens to be equal to lim/ZA(- | er), when Li,L 2 ,L 3 — > 
oo, under a fixed boundary condition a, we shall call it the Gibbs state 
associated to the boundary condition a. One says also that this state exists 
in the thermodynamic limit. Then, equivalently, the correlation functions 
(|4]) converge to the corresponding expectation values in this state. 

This model presents, at low temperatures, that is, for (3 > f3 c , where (3 C is 
the critical inverse temperature, two different thermodynamic pure phases, a 
dense and a dilute phase in the lattice gas language (called here the positive 
and the negative phase). This means two extremal translation invariant 
Gibbs states, /z + and obtained as the Gibbs states associated with the 
boundary conditions a, respectively equal to the ground configurations a(i) = 
1 and a(i) = — 1, for all i G C. The spontaneous magnetization, 



is then strictly positive. On the other hand, if /3 < (3 C , then the Gibbs state 
is unique and m* = 0. 

Each configuration inside A can be described in a geometric way by spec- 
ifying the set of Peierls contours which indicate the boundaries between the 
regions of spin 1 and the regions of spin —1. Unit square faces are placed 
midway between the pairs of nearest-neighbor sites i and j, perpendicularly 
to these bonds, whenever a(i)a(j) — — 1. The connected components of this 
set of faces are the Peierls contours. Under the boundary conditions (+) 
and (— ), the contours form a set of closed surfaces. They describe the de- 
fects of the considered configuration with respect to the ground states of the 




(4) 



H (a(i)), 



(5) 
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system (the constant configurations 1 and —1), and are a basic tool for the 
investigation of the model at low temperatures. 

In order to study the interface between the two pure phases one needs 
to construct a state describing the coexistence of these phases. This can be 
done by means of a new boundary condition. Let n = (ni,n 2 ,n 3 ) be a unit 
vector in R 3 , such that n 3 > 0, and introduce the mixed boundary condition 
(±, n), for which 

- / -\ f 1 if % ■ n > 0, lR , 
^H-l if i-n<0. (6) 
This boundary condition forces the system to produce a defect going trans- 
versally through the box A, a big Peierls contour that can be interpreted as 
the microscopic interface (also called a domain wall). The other defects that 
appear above and below the interface can be described by closed contours 
inside the pure phases. 

The free energy, per unit area, due to the presence of the interface, is the 
surface tension. It is defined by 

r(n)= hm lim _-^ln^^ (7) 

In this expression the volume contributions proportional to the free energy 
of the coexisting phases, as well as the boundary effects, cancel, and only the 
contributions to the free energy due to the interface are left. The existence 
of such a quantity indicates that the macroscopic interface, separating the 
regions occupied by the pure phases in a large volume A, has a microscopic 
thickness and can therefore be regarded as a surface in a thermodynamic 
approach. 

Theorem 1 The interfacial free energy per unit area, r(n) ; exists, is boun- 
ded, and its extension by positive homogeneity, /(x) = |x| r(x/|x|) 7 is a 
convex function on R 3 . Moreover, r(n) is strictly positive for (3 > f3 c , and 
vanishes if (3 < (3 C . 

The existence of r(n) and also the last statement were proved by Lebowitz 
and Pfister (1981), in the particular case n = (0,0,1), with the help of 
correlation inequalities. A complete proof of the theorem was given later 
with similar arguments. The convexity of / is equivalent to the fact that the 
surface tension r satisfies a thermodynamic stability condition known as the 
pyramidal inequality (see Messager et al, 1992). 
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3 Gibbs states and interfaces 

In this section we consider the (±, n ) boundary condition, also simply de- 
noted (±), associated to the vertical direction n = (0,0, 1), 

a(i) = 1 if i 3 > 0, a(i) = -1 if i 3 < 0. (8) 

The corresponding surface tension is = r(n ). We shall first recall some 
classical results which concern the Gibbs states and interfaces at low tem- 
peratures. 

According to the geometrical description of the configurations introduced 
in Section 2, we observe that 

Z^(A)/Z+(A) = ]Texp ( - 2/3J\X\ - U A (X)), (9) 

A 

where the sum runs over all microscopic interfaces A compatible with the 
boundary condition and |A| is the number of faces of A (inside A). The 
term U\(X) equals — In Z + (X, A)/Z + (A), the sum in the partition function 
Z + (X, A) being extended to all configurations whose associated contours do 
not intersect A. Each term in sum (Q gives a weight proportional to the 
probability of the corresponding microscopic interface. 

At low (positive) temperatures, we expect the microscopic interface cor- 
responding to this boundary condition, which at zero temperature coincides 
with the plane i 3 = —1/2, to be modified by small deformations. Each mi- 
croscopic interface A can then be described by its defects, with respect to 
the interface at (3 = oo. To this end one introduces some objects, called 
walls, which form the boundaries between the horizontal plane portions of 
the microscopic interface, also called the ceilings of the interface. 

More precisely, one says that a face of A is a ceiling face if it is horizontal 
and such that the vertical line passing through its center does not have other 
intersections with A. Otherwise, one says that it is a wall face. The set of wall 
faces splits into maximal connected components. The set of walls, associated 
to A, is the set of these components, each component being identified by its 
geometric form and its projection on the plane is = —1/2. Every wall u, with 
projection ir(u), increases the energy of the interface by a quantity 2J||w||, 
where ||w|| = \ou\ — \tt(uj)\, and two walls are compatible if their projections 
do not intersect. In this way the microscopic interfaces may be interpreted 
as a "gas of walls" on the two-dimensional lattice. 
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Dobrushin, who developed the above analysis, proved also the dilute char- 
acter of this "gas" at low temperatures. This implies that the microscopic 
interface is essentially flat, or rigid. One can understand this fact by noticing 
first, that the probability of a wall is less than exp(— 2(5 J||a;||), and second, 
that in order to create a ceiling in A, which is not in the plane i% = —1/2, 
one needs to surround it by a wall, that one has to grow when the ceiling is 
made over a larger area. 

Using correlation inequalities one proves that the Gibbs state /r* 3 , associ- 
ated to the (±) boundary conditions, always exists, and that it is invariant 
under horizontal translations of the lattice, i.e., u ± {a{A + a)) = n ± {a{A)) for 
all a = (ai,ci2,0). It is also an extremal Gibbs state. Let m(z) be the mag- 
netization /^((^(z)) at the site z = (0, 0, z). The function m(z) is monotone 
increasing and satisfies the symmetry property m(—z) = —m(z + 1). Some 
consequences of Dobrushin's work are the following properties. 

Theorem 2 If the temperature is low enough, i.e., if [5J > c\, where C\ is a 
given constant, then 



Equation ( FlOi) is just another way of saying that the interface is rigid 
and that the state u, is non-translation invariant (in the vertical direction). 
Then, the correlation functions /i ± (o"(A)) describe the local properties, or 
local structure, of the macroscopic interface. In particular, the function m(z) 
represents the magnetization profile. Then statement (Till) , together with the 
symmetry property, tells us that the thickness of this interface is finite, with 
respect to the unit lattice spacing. 

The statistics of interfaces has been rewritten in terms of a gas of walls 
and this system may further be studied by cluster expansion techniques. 
There is an interaction between the walls, coming from the term Ua(X) in 
equation (Q, but a convenient mathematical description of this interaction 
can be obtained by applying the standard low temperature cluster expansion, 
in terms of contours, to the regions above and below the interface. 

This method was introduced by Gallavotti in his study, mentioned below, 
of the two-dimensional Ising model. It has been applied by Bricmont et ai, 
to examine the interface structure in the present case. As a consequence, it 
follows that the surface tension, more exactly /3r ± (/3), and also the correlation 



m (0) is strictly positive, 

m ± (z) — > m*, when z — > oo, exponentially fast. 



(10) 
(11) 
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functions, are analytic functions at low temperature. They can be obtained 
as explicit convergent series in the variable £ = e~ 2/3J . 

The same analysis applied to the two-dimensional model shows a very 
different behavior at low temperatures. In this case, the microscopic interface 
A is a polygonal line and the walls belong to the one-dimensional lattice. One 
can then increase the size of a ceiling without modifying the walls attached 
to it. 

Indeed, Gallavotti turned this observation into a proof that the Gibbs 
state ^ is now translation invariant. The line A undergoes large fluctuations 
of order y/Li, and disappears from any finite region of the lattice, in the 
thermodynamic limit. In particular, we have then /x ± = (1/2) + a 
result that extends to all boundary conditions (±, n). 

Using these results Bricmont et al. studied also the local structure of 
the interface at low temperatures and showed that its intrinsic thickness is 
finite. To study the global fluctuations one can compute the magnetization 
profile by introducing, before taking the thermodynamic limit, a change of 
scale: ^(a(zL{)), with 5 = 1/2 or near to this value. This is an exact 
computation that has been done by Abraham and Reed. 

Let us come back to the three-dimensional Ising model where we know 
that the interface orthogonal to a lattice axis is rigid at low temperatures. 

Question 1 At higher temperatures, but before reaching the critical temper- 
ature, do the fluctuations of this interface become unbounded, in the thermo- 
dynamic limit, so that the corresponding Gibbs state is translation invariant? 

One says then that the interface is rough, and it is believed that, effec- 
tively, the interface becomes rough when the temperature is raised, undergo- 
ing a roughening transition at an inverse temperature /3r > /3 C . 

It is known that (3r < (5f =2 , the critical inverse temperature of the two- 
dimensional Ising model, since van Beijeren proved using correlation inequali- 
ties, that above this value, the state ^ is not translation invariant. Recalling 
that the rigid interface may be viewed as a two-dimensional system, the sys- 
tem of walls, a representation that would become inappropriate for a rough 
interface, one might think that the phase transition of the two-dimensional 
Ising model is relevant for the roughening transition, and that [5r is some- 
where near (3f =2 . Indeed, approximate methods, used by Weeks et al, give 
some evidence for the existence of such a (3r and suggest a value slightly 
smaller than (3^ =2 , as shown in Table 1. To this day, however, it appears to 
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be no proof of the fact that /3r > j3 c , i.e., that the roughening transition for 
the three-dimensional Ising model really occurs. 



Table 1. Some temperature values. 



d = 3 (3 C J ~ 0.22 approximate critical temperature 
d = 3 (3 R J ~ 0.41 conjectured roughening temperature 
d = 2 (3 C J = 0.44 exact critical temperature 



At present one is able to study rigorously the roughening transition only 
for some simplified models with a restricted set of admissible microscopic 
interfaces. Moreover, the closed contours, describing the defects above and 
below A, are neglected, so that these two regions have the constant configu- 
rations 1 or —1, and one has U\(X) = in equation (J5J). 

The best known of these models is the classic SOS (solid-on-solid) model 
in which the interfaces A have the property of being cut only once by all 
vertical lines of the lattice. This means that A is the graph of a function 
that can equivalently be used to define the possible configurations of A. If 
A contains the horizontal face with center (ii,i2,h — 1/2), then the value at 
(hjiz) of the associated function is <j>(ii,i2) = h- 

The proof that the SOS model with the boundary condition (±) has 
a roughening transition is a highly non-trivial result due to Frohlich and 
Spencer. When is small enough, the fluctuations of A are of order Vln L 
(in a cubic box of side L). 

Moreover, other interface models, with additional conditions on the al- 
lowed microscopic interfaces, are exactly solvable. The BCSOS (body-cen- 
tered solid-on-solid) model, introduced by van Beijeren, belongs to this class. 
It is, in fact, the first model for which the existence of a roughening tran- 
sition has been proved. More recently, also the TISOS (triangular Ising 
solid-on-solid) model, introduced by Blote and Hilhorst and further studied 
by Nienhuis et at, has been considered in this context. 

The interested reader can find more information and references, concern- 
ing the subject of this section, in the review article by Abraham (1986). 



10 



S. Miracle-Sole 



4 Wetting phenomena 

Next we consider the Ising model over a plane horizontal substrate (also 
called a wall) and study the difference of surface tensions which governs the 
wetting properties of this system. 

We first describe the approach developed by Frohlich and Pfister (1987) 
and briefly report some results of their study. We consider the model on the 
semi-infinite lattice 

C' = {i e Z 3 : i 3 > 0}. (12) 

A magnetic field, K > 0, is added on the boundary sites, 23 = 0, which 
describes the interaction with the substrate, supposed to occupy the comple- 
mentary region C\C. 

We constrain the model in the finite box A' = A n C, with A as above, 
and impose the value of the spins outside. The hamiltonian becomes 

#A'(°A' \&) = —J £ a{i)a{j)-K £ o(i). (13) 
(tj)nAV0 ieA',« 3 =o 

Here <ta' represents the configuration inside A', the pairs are contained 
in £', and a(i) = a(i) when % (jL A', the configuration a being the given 
boundary condition. The corresponding partition function is denoted by 
Z W °(A'). 

Since there are two pure phases in the model we must consider two surface 
free energies, or surfaces tensions, r w+ and r w_ , between the wall and the 
positive or negative phase present in the bulk. They are defined through the 
choice of the boundary condition, a(i) = 1 or a(i) = —1, for all i G C. Let 
us consider first the case of the (— ) boundary condition. 

The surface free energy contribution, per unit area, due to the presence 
of the wall, when we have the negative phase in the bulk, is 

T w -(fi,K)= lim lim In K ~4 (14) 

The division by Z~(A)5 allows us to subtract from the total free energy, 
lnZ w ~(A'), the bulk term and all boundary terms which are not related to 
the presence of the wall. The existence of limit (fT4"]) follows from correlation 
inequalities, and we have r w ~ > 0. 
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Figure 2. Boundary conditions for the cubic lattice. Above, the box A with 
the (±) and (step) boundary conditions. Below, the box A' and the wall W 
with the (w— ) boundary conditions. 



One can prove, as well, the existence of the Gibbs state fi w ~ of the semi- 
infinite system, associated to the (— ) boundary condition. This state is the 
limit of the finite volume Gibbs measures /iA'(°A' | (— )) defined by hamil- 
tonian ffT3]) . It describes the local equilibrium properties of the system near 
the wall, when deep inside the bulk the system is in the negative phase. 
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Similar definitions give the surface tension r w+ and the Gibbs state /i w+ , 
corresponding to the boundary condition a(i) = 1, for all i G A'. 

We remark that the states fi w+ and /i w ~ are invariant by translations 
parallel to the plane 13 = 0, and introduce the magnetizations, m w ~(z) = 
fi w ~(<j(z)), where z denotes the site (0,0, z), m w_ = m w ~(0), and similarly 
rn f!+ {z) and m w+ . Their connection with the surface free energies is given 
by the formula 

t w -{P,K) -t w+ (P,K) = I (m w+ (/3,s) -m w -([3,s))ds. (15) 

J 

We mention in the following theorem some results of Frohlich and Pfister's 
study. Here r ± is, as before, the usual surface tension between the two pure 
phases of the system, for a horizontal interface. 

Theorem 3 With the above definitions, we have 

t-(P,K)-t- + (P,K) <T±{p), (16) 
m w+ (/3, K) - m w -(/3, K) > 0, (17) 



and the difference in (11) is a monotone decreasing function of the parameter 



K. Moreover, if m w+ = m w , then the Gibbs states fi w+ and /i w coincide. 

The proof is a subtle application of correlation inequalities. Since, from 
Theorem [31 the integrand in equation (fT5|) is a positive and decreasing func- 
tion, the difference At = r w ~ — r w+ is a monotone increasing and concave 
(and hence continuous) function of the parameter K. On the other hand, 
one can prove that At = if K > J. This justifies the following definition 

KM = min{K: At((3,K) = r + ~ (/?)}. (18) 

In the thermodynamic description of wetting, the partial wetting regime 
is characterized by the strict inequality in equation ffTB]) . Equivalently, by 
K < K w (f3). We must have then m w+ 7^ m w_ , because of equation ([151) . 
This shows that, in the case of partial wetting, /z w+ and /i w ~ are different 
Gibbs states. 

The complete wetting regime is characterized by the equality in equation 
fTTol) . that is, by K > K w (f3). Then, we have m w+ = m w ~, and taking 
into account the last statement in Theorem [31 also /i w+ = fi w ~~ . This last 
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result implies that there is only one Gibbs state. Thus complete wetting 
corresponds to unicity of the Gibbs state. 

In this case, we also have lim m w ~(z) = m*, when z — >■ oo, because this 
is always true for m w+ (z). This indicates that we are in the positive phase 
of the system although we have used the (— ) boundary condition, so that 
the bulk negative phase cannot reach anymore the wall. The film of positive 
phase, which wets the wall completely, has an infinite thickness with respect 
to the unit lattice spacing, in the thermodynamic limit. 

When (3 = oo only few particular ground configurations contribute to the 
partition functions, such as the configuration a{i) = —1 for the partition 
function Z w ~, etc., and we obtain At = 2K and = 2 J. For non zero 
but low temperatures the small perturbations of these ground states have 
to be considered, a problem that can be treated by the method of cluster 
expansions. In fact, the corresponding defects can be described by closed 
contours as in the case of pure phases. 

Theorem 4 For K < J, the functions (3t w ~((3,K) and /3t w+ ((3,K) are 
analytic at low temperatures, i.e., provided that (3 ( J — K) > C2, where Ci is 
a given constant. Moreover, rri w+ (z) and m w_ (z) tend, respectively, to m* 
and to —m*, when z — > oo, exponentially fast. 

The last statement in Theorem H] tells us that the wall affects only a layer 
of finite thickness (with respect to the lattice spacing). From a macroscopic 
point of view the negative phase reaches the wall, and we are in the partial 
wetting regime. Indeed, a strict inequality holds in equation (JTB)) . 

Thus, for K < J there is always partial wetting at low temperatures. 
Then the following question arises: 

Question 2 Is there a situation of complete wetting at higher temperatures ? 
It is understood here that K takes a fixed value, characteristic of the substrate, 
such that < K < J . 

This is known to be the case in dimension d = 2, where the exact value 
of K w ((3) can be obtained from Abraham's solution of the model: 

cosh 2f3K w = cosh 2/3 J - e" 2/3J sinh 2(3 J. 

Then complete wetting occurs for (3 in the interval (3 C < (3 < f3 w (K), where (3 C 
is the critical inverse temperature and (3 W (K) is the solution of K w {(3) = K. 
The case d = 2 has been reviewed in Abraham (1986). 
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To our knowledge, the above question remains an open problem for the 
Ising model in dimension d = 3. The problem has however been solved for 
the simpler case of a SOS interface model. In this case, a nice and rather 
brief proof has been given by Chalker (1982) of the following result: One has 
m w+ = m w_ , and hence complete wetting, if 

2/3(J- K) < -ln(l -e" 8/3J ). 

It is very plausible that a similar statement is valid for the semi-infinite 
Ising model and, also, that Chalker's method could play a role for extending 
the proof to this case, provided an additional assumption is made. Namely, 
that is sufficiently large, and hence J — K small enough, in order to insure 
the convergence of the cluster expansions and to be able to use them. 

5 Equilibrium crystals 

The shape of an equilibrium crystal is obtained, according to thermody- 
namics, by minimizing the surface free energy between the crystal and the 
medium, for a fixed volume of the crystal phase. Given the orientation de- 
pendent surface tension r(n), the solution to this variational problem, known 
under the name of Wulff construction, is the following set: 

W = {x G R 3 : x • n < r(n) for all n}. (19) 

Notice that the problem is scale invariant, so that if we solve it for a given 
volume of the crystal, we get the solution for other volumes by an appropriate 
scaling. We notice also that the symmetry r(n) = r(— n) is not required for 
the validity of formula ( TT9l) . In the present case r(n) is obviously a symmetric 
function, but non-symmetric situations are also physically interesting and 
appear, for instance, in the case of a drop on a wall discussed in Section 4. 

The surface tension in the Ising model, between the positive and negative 
phases, has been defined in equation ([7j). In the two-dimensional case, this 
function r(n) has (as shown by Abraham) an exact expression in terms of 
some Onsager's function. It follows (as explained in Miracle-Sole, 1999) that 
the Wulff shape W, in the plane (x\, x 2 ), is give by 

cosh/3xi + cosh/3x 2 < cosh 2 2/3 J / sinh2/3J. 

This shape reduces to the empty set for (3 < (3 C , since the critical (3 C satisfies 
sinh2J/3 c = 1. For > (3 C , it is a strictly convex set with smooth boundary. 
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Figure 3. Cubic equilibrium crystal shown in a projection parallel to the 
(1,1,1) direction. The three regions 1), 2) and 3) indicate the facets and the 
remaining area represents a curved part of the crystal surface. 

In the three-dimensional case, only certain interface models can be ex- 
actly solved (see Section 3). Consider the Ising model at zero temperature. 
The ground configurations have only one defect, the microscopic interface A, 
imposed by the boundary condition (±, n). Then, from equation flHJ), we may 
write 

r(n) = lim (£ A (n) - /T^n)), (20) 

where E\ = 2J\X\ is the energy (all A have the same minimal area) and 
]Va the number of the ground states. Every such A has the property of 
being cut only once by all straight lines orthogonal to the diagonal plane 
h + ^2 + H = 0, provided that n k > 0, for k = 1,2, 3. Each A can then be 
described by an integer function defined on a triangular plane lattice, the 
projection of the cubic lattice £ on the diagonal plane. The model defined 
by this set of admissible microscopic interfaces is precisely the TISOS model. 
A similar definition can be given for the BCSOS model that describes the 
ground configurations on the body-centered cubic lattice. 

From a macroscopic point of view, the roughness or the rigidity of an 
interface should be apparent when considering the shape of the equilibrium 
crystal associated with the system. A typical equilibrium crystal at low 
temperatures has smooth plane facets linked by rounded edges and corners. 
The area of a particular facet decreases as the temperature is raised and 
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the facet finally disappears at a temperature characteristic of its orientation. 
It can be argued that the disappearance of the facet corresponds to the 
roughening transition of the interface whose orientation is the same as that 
of the considered facet. 

The exactly solvable interface models mentioned above, for which the 
function r(n) has been computed, are interesting examples of this behav- 
ior, and provide a valuable information on several aspects of the roughening 
transition. This subject has been reviewed by Abraham (1986), van Beijeren 
and Nolden (1987), and Kotecky (1989). 

For example, we show in figure 2 the shape predicted by the TISOS 
model (one eighth of the shape because of the condition n& > 0). In this 
model the interfaces orthogonal to the three coordinate axes are rigid at low 
temperatures. 

For the three-dimensional Ising model at positive temperatures, the de- 
scription of the microscopic interface, for any orientation n, appears as a 
very difficult problem. It has been possible, however, to analyze the in- 
terfaces which are very near to the particular orientations n , discussed in 
Section 3. This analysis allows us to determine the shape of the facets in a 
rigorous way. 

We first observe that the appearance of a facet in the equilibrium crystal 
shape is related, according to the Wulff construction, to the existence of 
a discontinuity in the derivative of the surface tension with respect to the 
orientation. More precisely, assume that the surface tension satisfies the 
convexity condition of Theorem [TJ and let this function r(n) = t(9, 0) be 
expressed in terms of the spherical coordinates of n, the vector n being 
taken as the x 3 axis. A facet orthogonal to n appears in the Wulff shape if, 
and only if, the derivative dr(9,(f))/d9 is discontinuous at the point 9 = 0, 
for all <p. The facet J 7 C dW consists of the points x e R 3 belonging to the 
plane x^ = r(n ) and such that, for all <j) between and 2n, 

Xi cos (j) + x 2 sin (j) < dr(9, 4>)/d9\ e=0 + (21) 

The step free energy is expected to play an important role in the facet 
formation. It is defined as the free energy associated with the introduction of 
a step of height 1 on the interface, and can be regarded as an order parameter 
for the roughening transition. Let A be a parallelepiped as in Section 2, and 
introduce the (step, m) boundary conditions, associated to the unit vectors 
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m = (cos 0, sin 0) 6 R , by 

_/.x_ fl if i>0 or if i 3 = and i±mi + i 2 m2 > 0, 
I —1 otherwise. 

Then, the step free energy, per unit length, for a step orthogonal to m (with 
m 2 > 0) on the horizontal interface, is 

r step (0) = lim lim lim In ^'"a^ (23) 

A first result concerning this point, was obtained by Bricmont et ai, by 
proving a correlation inequality which establish r stcp (0) as a lower bound to 
the one-sided derivative dr(9,0)/d9 at 9 = + (the inequality extends also 
to <fi 7^ 0). Thus when r stcp > a facet is expected. 

Using the perturbation theory of the horizontal interface, it is possible to 
study also the microscopic interfaces associated with the (step, m) boundary 
conditions. When considering these configurations, the step may be viewed 
as an additional defect on the rigid interface described in Section 2. It is, in 
fact, a long wall going from one side to the other side of the box A. The step 
structure at low temperatures can then be analyzed with the help of a new 
cluster expansion. As a consequence of this analysis we have the following 
theorem. 

Theorem 5 If the temperature is low enough, i.e., if f3J > c 3; where c 3 is a 
given constant, then the step free energy, r step (0), exists, is strictly positive, 
and extends by positive homogeneity to a strictly convex function. More- 
over, /3r stcp (0) is an analytic function of ( = e _2J/3 , for which an explicit 
convergent series expansion can be found. 

Using the above results on the step structure, similar methods allow us 
to evaluate the increment in surface tension of an interface titled by a very 
small angle 9 with respect to the rigid horizontal interface. This increment 
can be expressed in terms of the step free energy and one obtains the following 
relation. 

Theorem 6 For f3J > c^, we have 

dr(M)/^U =o+ =r step (0). (24) 
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This relation, together with equation (121j) . implies that one obtains the 
shape of the facet by means of the two-dimensional Wulff construction applied 
to the step free energy. The reader will find a detailed discussion on these 
points, as well as the proofs of Theorems [5] and O in Miracle-Sole (1995). 

From the properties of r stcp stated in Theorem [5l it follows that the 
Wulff equilibrium crystal presents well defined boundary lines, smooth and 
without straight segments, between a rounded part of the crystal surface and 
the facets parallel to the three main lattice planes. 

It is expected, but not proved, that at a higher temperature, but be- 
fore reaching the critical temperature, the facets associated with the Ising 
model undergo a roughening transition. It is then natural to believe that the 
equality (1241 is true for any (3 larger than (3r, allowing us to determine the 
facet shape from equations (I21~j) and (|24|) . and that for (3 < (3r, both sides 
in this equality vanish, and thus, the disappearance of the facet is involved. 
However, the condition that the temperature is low enough is needed in the 
proofs of Theorems and El 

See also: Gibbs states. Ising model. Phase transitions. Wulff droplets. 
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